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SUMMARY 


/ 


The  notion  of  the  driving  traction  on  a  surface  of  strain  discontinuity 
in  a  continuum  undergoing  a  general  thermo -mechanical  process  is  defined 
and  discussed.  In  addition,  the  associated  constitutive  notion  of  a  kinetic 
XeJL&tjpfl,  in  which  the  normal  velocity  of  propagation  of  the  surface  of 
discontinuity  may  be  a  given  function  of  the  driving  traction  and  tempera¬ 
ture,  is  introduced  for  the  special  case  of  a  thermoelastic  material . 

V 
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1.  Introduction 

Various  aspects  of  the  theory  of  finite  elastostatics  for  materials 
characterized  by  non- elliptic  elastic  potentials  have  been  studied  in  a 
number  of  recent  papers# ^ee  for  example  Abeyaratne  and  Knowles  (1987a, b, 
1988a, b),  Ball  (1977)  ,  Bali  and  James  (1987),  Ericksen  (1975),  Fosdick  and 
MacSithigh  (1983),  Gurtin  (1983),  James  (1981,1986),  Knowles  and  Sternberg 
(1978),  Silling  (1988).  One  feature  of  such  materials  is  that,  under  suit¬ 
able  conditions,  they  can  sustain  equilibrium  deformations  in  which  the 
displacement  gradient  and  stress  tensors  suffer  jump  discontinuities  across 
certain  surfaces  in  the  body,^  while  the  displacement  and  traction  remain 
continuous;  such  singular  surfaces  have  been  called  equilibrium  shocks.  One 
area  in  which  the  associated  theory  finds  application  is  that  of  the  conti¬ 
nuum-mechanical  modelin  g/of  a  solid  in  equilibrium  with  more  than  one 
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"phase"  present.  In  ptls  setting,  an  equilibrium  shock  corresponds  to  the 
boundary  between  two  distinct  phases  of  the  material  (Ball  and  James 
(1987),  Janwrs  (1981,1986),  Silling  (1988)). 


A  quasi-static  motion  involving  equilibrium  shocks  at  each  instant  may 
be  dissipative  (Knowles  (1979)):  in  any  portion  of  the  body  that  is  trav¬ 
ersed  by  a  moving  shock,  the  rate  of  work  of  the  external  forces  differs 
from  the  rate  of  storage  of  strain  energy  by  the  rate  of  work  done  in  mov¬ 
ing  the  surface  of  discontinuity.  This  latter  rate  of  work  can  be  expressed 
as  the  integral  over  the  shock  surface  of  the  product  of  a  scalar  "driving 
traction"  f  with  the  component  of  shock  velocity  normal  to  the  shock 
itself.  At  each  instant  during  the  motion,  the  value  of  f  at  each  point  on 
the  singular  surface  can  be  calculated  from  the  limiting  values  of  the 
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deformation  gradient  on  the  two  sides  of  the  surface.  If  the  driving  trac¬ 
tion  f  vanishes  at  all  points  of  the  shock  for  each  equilibrium  state 
visited  during  a  quasi-static  motion,  the  motion  is  dissipation- free .  An 
equilibrium  state  for  which  f  -  0  is  said  to  satisfy  the  Maxwell  condition. 

In  the  presence  of  equilibrium  shocks,  the  mechanical  balance  laws  lead 
to  the  usual  field  equations  at  points  in  the  body  away  from  the  singular 
surfaces  and  require  the  traction  to  be  continuous  across  these  surfaces. 
However,  on  solving  even  the  simplest  of  boundary -value  problems  for  non- 
elliptic  elastic  materials,  one  encounters  a  massive  failure  of  uniqueness 
of  solution  (Abeyaratne  (1980),  Abeyaratne  and  Knowles  (1987a)).  One  way  to 
single  out  a  preferred  equilibrium  field  from  among  the  infinitely  many 
available  ones  is  to  require  that  the  field  be  stable  in  the  sense  that  the 
associated  potential  energy  be  an  absolute  minimum.  If  an  equilibrium  field 
containing  a  shock  is  to  be  stable  in  this  sense,  it  is  known  that  the  Max¬ 
well  condition  f  -  0  Is  necessary  (Abeyaratne  (1983),  Ericksen  (1975),  Gur- 
tin  (1983),  James  (1981)). 

A  somewhat  more  general  point  of  view  holds  that  the  lack  of  uniqueness 
in  the  conventional  equilibrium  problems  for  non-elliptic  elastic  materi¬ 
als  arises  from  a  constitutive  deficiency  associated  with  particles  on  the 
shock  surface.  Such  a  view  was  adopted  by  Abeyaratne  and  Knowles  (1988a, b) 
in  the  one -dimensional  context  of  bar  theory.  By  exploiting  an  analogy 
between  the  problem  considered  and  internal -variable  theories  of  inelastic 
solids  (see,  for  example,  Rice  (1971,1975)),  they  are  led  to  postulate  a 
supplementary  constitutive  requirement  in  the  form  of  a  "kinetic  relation" 
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be tween  f  and  the  velocity  of  the  shock  during  a  quasi-static  motion.  This 
leads  to  a  determinate  macroscopic  response  (or  force-elongation  relation) 
for  the  bar  in  quasi-static  motions.  In  general,  this  response  exhibits 
rate-  and  his tory- dependence .  Two  limiting  cases  of  the  kinetic  relation 
describe  rate -independent  behavior:  one  corresponds  to  the  Maxwell  condi¬ 
tion  and  hence  to  dissipation-free,  reversible  macroscopic  response.  The 
other  leads  to  a  force-elongation  relation  similar  to  that  associated  with 
rate - independent  plasticity. 

The  investigations  described  above  are  for  the  most  part  carried  out 
within  the  framework  of  the  purely  mechanical  theory  of  non-elliptic  elas¬ 
tic  materials.  They  are  limited  to  the  study  of  equilibrium  states  or 
one-parameter  families  of  such  states  (i.e.,  quasi-static  motions),  and 
they  are  often  confined  to  one -dimensional  settings  as  well.  The  purpose  of 
the  present  study  is  to  consider  the  corresponding  issues  in  a  more  gen¬ 
eral,  three-dimensional  setting  in  which  the  material  need  not  be  elastic, 
and  both  thermal  and  inertial  effects  are  taken  into  account.  In  Section  2, 
we  recall  the  basic  thermodynamic  and  mechanical  laws  and  the  associated 
field  equations,  inequalities  and  jump  conditions  for  a  continuum  under¬ 
going  a  thermo -mechanical  process  in  which  displacement  and  temperature  are 
assumed  to  be  continuous,  but  their  gradients  are  permitted  to  jump  across 
a  moving  surface.  We  derive  a  useful  representation  of  the  entropy  produc¬ 
tion  rate  in  Section  3,  and  we  use  it  to  introduce  the  notion  of  a  driving 
traction  acting  on  a  singular  surface  in  an  arbitrary  continuum  during  a 
thermo -mechanical  process  of  the  assumed  type.  We  specialize  the  earlier 
results  to  isothermal  dynamic  processes  in  Section  4.  No  constitutive 
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assumptions  are  invoked  until  Section  5,  where  we  specialize  the  foregoing 
results  to  the  case  of  a  continuum  composed  of  a  thermo -elastic  material. 
Section  6  is  devoted  to  a  discussion  of  the  notion  of  a  constitutive 
“kinetic  relation" . 

2.  Balance  laws,  field  equations  and  jump  conditions 

Consider  a  body  B  that  occupies  a  region  R  in  a  reference  configura¬ 
tion.  A  motion  of  the  body  on  a  time  interval  [tQ.tjJ  is  characterized  by 
a  one-parameter  family  of  invertible  mappings  y(-,t):  R  -*■  R*  ,  with 

y  -  y(x,t)  -  x  +  u(x, t)  for  X6R,  te[t0,ti].  (2.1) 

We  assume  that  the  deformation  y,  or  equivalently  the  displacement  u_,  is 
continuous  with  piecewise  continuous  first  and  second  derivatives  on 
RxftQ.t^].  Let  F(x ,  t)  -  Grad£(x,  t)  and  v(x, t)  -  d£(£,t)/3t  stand 
respectively  for  the  deformation  gradient  tensor  and  the  particle  velocity 
at  points  (x,t)  in  space-time  where  they  exist. 


Let  p(x)  denote  the  mass  density  of  B  at  the  point  x  in  the  reference 
configuration,  b(x,t)  the  body  force  per  unit  mass,  and  a(x,t)  the  nominal 
stress  tensor.  At  each  t,  we  require  p(-)  and  b(-,t)  to  be  continuous  on  R, 
while  j7(-,t)  is  to  be  piecewise  continuous  with  a  piecewise  continuous 
gradient  on  R.  The  balance  laws  for  linear  and  angular  momentum  require 


• 

» 

cm  dA  + 

pb  dV  -  d/dt 

3D 

D 

pv  dV  , 


(2.2) 
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0  x  an  dLA  + 

x  «•> 

3D 


x  pb  dV 


d/dt 


A 

y  x  pv  dV 


(2.3) 


respectively,  at  each  te[tQ,t^]  and  for  all  regular  subregions  DcR. 

Next,  let  q(x,t)  denote  the  nominal  heat  flux  vector,  r(x,t)  the  heat 
supply  per  unit  mass  and  c(x,t)  the  internal  energy  per  unit  mass.  At  each 
t,  we  suppose  that  r(-,t)  is  continuous  on  R  and  that  q(-,t)  is  piecewise 
continuous  with  a  piecewise  continuous  gradient  on  R.  The  internal  energy 
c (•,•)  is  required  to  be  piecewise  continuous  with  piecewise  continuous 
first  derivatives  on  Rx[tg,t^].  The  first  law  of  thermodynamics  requires 
that  at  each  instant  t, 


• 

• 

* 

on*v  dA  + 

^b*v  dV  + 

q*n  dA  + 

3D 

• 

D 

3D 

d/dt 


pr  dV 


pc  dV  +  d/dt 


(1/2)  pv.v  dV,  (2.4) 

D 


for  every  regular  DCR.  Finally,  let  0(x,t)  denote  the  absolute  temperature 
and  »?(x,t)  the  entropy  per  unit  mass.  At  each  t,  we  assume  that  0(-,t)  is 
continuous  with  a  piecewise  continuous  gradient  on  R,  while  »?(•,•)  is 
assumed  to  be  piecewise  continuous  with  piecewise  continuous  first  deriva¬ 
tives  on  Rx[tg,t]J.  The  rate  of  entropy  production  in  a  regular  region  DcR 
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is  defined  to  be 


r(t;D)  -  d/dt  |  pt)  dV 
D 


q*n /9  dA  -  j  pr/0  dV. 
3D  D 


(2.5) 


The  Clausius -Duhem  version  of  the  second  law  of  thermodynamics  requires 


r(t;D)  >  0  for  DcR,  te[t0,t;|J. 


(2.6) 


At  a  fixed  instant  t,  localization  of  the  balance  laws  (2.2)-(2.4)  and 

the  inequality  (2.6)  at  a  point  x  at  which  F,  F,  v,  q,  a,  c  and  rj  are  all 

^  ^ 

continuous  yields  the  following  familiar  local  results: 


Div  a  +  pb  -  pv, 


o*F  +  Div  q  +  pr  -  pc, 
Div(q/0)  +  pr/8  i  pfi . 


(2.7) 


On  the  other  hand,  suppose  that  S(t)  is  a  regular  surface  in  R  at  time  t 
across  which  some  or  all  of  the  thermo-mechanical  quantities  listed  above 
suffer  jump  discontinuities.  Localization  of  (2.2)-(2.6)  at  a  point  x.  on 
S(t)  yields  the  following  jump  conditions: 


[  Ign) )  +  P  [ [v] ]  Vn  -  0, 

[  [™»v]  ]  +  [[p(t  +  v*v/2)]]  Vn  +  [[q.n]]  -  0,  • 

tt Pi]]  Vn  +  [  [q-rj/^I  ]  <  0, 


on  S(t) , 


(2.8) 


where 


Vn  -  V-n,  (2.9) 

and  V  -  V(x,t)  is  the  velocity  of  the  point  x  on  the  moving  surface  S(t). 
The  unit  normal  n  on  the  singular  surface  S(t)  is  chosen  such  that  Vn  >  0; 
if  Vn  >  0,  the  positive  side  of  S(t)  is  the  side  into  which  (and  there¬ 
fore  n)  points.  If  g(x, t)  denotes  a  generic  field  quantity  that  jumps 
across  S(t),  we  write  [[g(x,  t)  ]  ]  -  g*(x,t)  -  g(x,t),  where  g(x,t)  and  g(x,t) 
stand  for  the  limiting  values  of  g  at  the  point  ^  on  S(t)  from  the  positive 
and  negative  sides,  respectively. 

In  addition  to  the  jump  conditions  listed  in  (2.8),  one  also  has  the 
kinematic  results 

[[FI]  1-  0  ,  [[v]]  -  -  [  [ F]  1  V  on  S(t)  ,  (2.10) 

where  is  any  vector  tangent  to  the  singular  surface  S(t).  The  jump 
conditions  in  (2.10)  are  immediate  consequences  of  the  smoothness  require¬ 
ments  imposed  on  the  deformation  (2.1). 

Conversely,  the  field  equations  (2.7)  together  with  the  jump  conditions 
(2.8),  (2.10)  imply  the  global  balance  laws  (2.2),  (2.3),  (2.4),  (2.6).  All 
of  the  above  results  may  be  found  in  Trusdell  and  Noll  (1965). 


In  the  present  section,  we  first  introduce  the  notion  of  the  shock 
driving  traction  associated  with  a  thermo -mechanical  process  in  an  arbi¬ 
trary  continuum.  The  process  is  assumed  to  possess  the  smoothness  specified 
in  the  preceding  section;  in  particular,  both  displacement  and  temperature 
are  required  to  be  continuous  across  the  moving  singular  surface  S(t).  We 
then  make  use  of  the  shock  driving  traction  to  derive  an  alternate  form  for 
the  entropy  and  energy  jump  conditions  (2.8)3  and  (2.8)2,  respectively. 

The  concept  of  shock  driving  traction  emerges  naturally  from  an  alter¬ 
nate  representation  for  the  rate  of  entropy  production  defined  in  (2.5).  To 
obtain  this  representation,  we  begin  by  putting  the  energy  jump  condition 
(2.8)2  into  a  form  more  useful  for  our  purposes.  In  the  identity 

[(an.v]]  -  (1/2)  (v  +  v) • [ [on]]  +  (1/2)  (tn  +  ojj)  •  [  [v]  J  ,  (3.1) 

we  replace  [ [on] ]  and  [ [ v] ]  on  the  right  by  substituting  from  (2.8)^  and 
(2.10)2,  respectively.  This  gives 

[[on-v]]  -  -(1/2)  p  [ [v»v] ]  Vn  -  (1/2)  (Jn  +£n).([[F]]  V)  .  (3.2) 

By  making  use  of  (2.10)^,  one  can  show  that 
(on  +  on) • ( [ [F] ]  V)  -  Vn  (|  +  o) • [ [F] ]  , 


(3.3) 


where  the  dot  product  on  the  right  is  that  associated  with  a  pair  of  ten¬ 
sors:  A*jB  -  Trace  (AB J) .  Combining  (3.2)  and  (3.3)  yields 

[Un-v]]  -  -  (1/2)  p  Vn  [  [v»v]  ]  -  (1/2)  Vn  (£  +  a).[[F]]  on  S(t).  (3.4) 

Using  (3.4)  to  replace  [ [gn*v] ]  in  (2.8)2  then  supplies  the  alternate 
version  of  the  energy  jump  condition: 

p  Vn  [Ml-  (1/2)  vn  (o  +  £)•[[£]]  -  [[q-n]]  onS(t).  (3.5) 


We  turn  now  to  (2.5)  and  assume  that,  at  each  instant  t,  the  region  DcR 
is  divided  into  two  parts  D+(t)  and  D'(t)  by  the  moving  surface  of  discon¬ 
tinuity  S(t).  Allowing  for  the  jump  in  specific  entropy  across  S(t),  we  may 
write  (2.5)  as 


k 

• 

» 

r(t;D)  - 

pr,  d V  - 

[[pr,]}  Vn  dA  - 

q»n /8  dA  - 

D  Si 

;t)nD 

3D  J 

pr/0  dV  , 


(3.6) 


which  is  equivalent  to 


r(t;D)  - 


pr\  d V  - 


[  [P*?]  )  Vn  dA  - 


q*n/0  dA  - 


S(t)oD 
px/8  dV  . 


- 

9DU3D 


[[q*n/5]]  dA 


S (t)OD 


(3.7) 
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Applying  the  divergence  theorem  to  the  third  integral  on  the  right  in 
(3.7)  and  then  utilizing  (2.7>3  to  eliminate  Div  q  yields 


r(t;D)  - 


• 

pv0  +  &  •  _F  -  pi 

+  q  •  Grad  8 

9 

82 

dV 


• 

4 

i  *  21 

]  +  [[pri9]]  Vn  ' 

8 

dA 


0.8) 


S(t)nD 


where  we  have  made  use  of  the  continuity  of  the  temperature  8  across 
S(t).  Finally,  version  (3.5)  of  the  energy  jump  condition  may  be  used  to 
eliminate  the  term  [[q«n]]  in  (3.8),  yielding  the  desired  representation 
for  the  rate  of  entropy  production: 

T(t;D)  -  rloc(t;D)  +  rcon(t;D)  +  rs(t;D)  ,  (3.9) 

where 


rioC<t;D) 


• .  •  • 

pr)8  +  a  •  F  -  pc 


9 


(3.10) 


rcon(t:D)  " 


q  •  Grad  8 


dV 


(3.11) 


ra(t;D)  - 


f  [  [  pc  -  pr)6  -  (£  +  a)  •  F/2  ]  ]  Vn  I 


dA  . 


(3.12) 


S(t)nD 


In  (3.9),  Che  CoCal  rate  of  entropy  production  T(t;D)  at  the  instant  t  for 
the  subregion  DcR  is  decomposed  into  three  parts:  r^oc  arises  from  local 
dissipation  in  the  material  away  from  the  singular  surface;  rcon  is  the 
entropy  production  rate  due  to  heat  conduction;  finally,  Ts  represents  the 
contribution  to  the  entropy  production  rate  arising  from  the  moving  singu¬ 
lar  surface  S(t).  A  similar  decomposition  in  the  absence  of  a  surface  of 
discontinuity  is  given  by  Truesdell  and  Noll  (1965),  §  79. 


The  internal  dissipation  6(x, t)  is  given  by 

6  m  dr)  +  (1  /p)a  *F  -  e  ;  (3.13) 

see  Chapter  2  of  Truesdell  (1969).  The  local  entropy  production  rate 
rioc(t;D)  of  (3.10)  may  be  written  in  terms  of  6  as  follows: 


I'loc^  ’ 


p(j)5(x,t)/« (x,t)  dV  . 
D 


(3.14) 


Next,  it  is  convenient  to  introduce  the  Helmholtz  free  energy  per  unit 
mass  rj>(i c,  t)  defined  by 

^(&,t)  -  c(x,t)  -  0(x,t)Tj(x,t)  ,  xeR,  teJtQ.t^]  .  (3.15) 
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We  now  define  the  (scalar)  shock  driving  traction  f(x, t)  on  the  singular 
surface  S(t)  by 

r  ' 

f(&,t)  -  p(jc)[[tf(£,t)]]  -(1/2)  g(x,t)+a(x,t)  *[[F(x,t)]]. 

k  4 

x€S(t),  teftQ.tjJ  .  (3.16) 

By  (3.15),  (3.16)  and  (3.12),  we  may  rewrite  the  contribution  Ts(t;D)  to 
the  rate  of  entropy  production  due  to  the  moving  singular  surface  in  terms 
of  the  shock  driving  traction  f(x, t),  the  temperature  0(x, t)  on  S(t)  and 
the  normal  velocity  Vn(jc,  t)  of  S(t): 

rs(t;D)  -  f(x,t)  Vn(x,t)/tf(x,t)  dA.  (3.17) 

S(t)nD 

We  note  from  (3.16)  that,  if  the  thermo -mechanical  process  under  con¬ 
sideration  is  smooth  in  the  sense  that  the  free  energy  ^  and  the  deforma¬ 
tion  gradient  are  continuous  everywhere  at  all  times ,  then  the  shock 
traction  f  on  any  surface  vanishes . 

Employing  the  representation  for  the  entropy  production  rate  furnished 
by  (3.9),  (3.11),  (3.14)  and  (3.17)  in  the  Clausius -Duhem  inequality  (2.6) 
and  localizing  the  result  at  a  point  x^ away  from  the  singular  surface  S(t) 
yields  the  inequality 

pS  +  (l/0)q*Grad  9  >  0  on  R-S(t),  te[t0>t1>].  (3.18) 
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This  result,  which  can  also  be  derived  directly  from  (2.7)3,  (2.7)^  and 
(3.13)  and  which  in  fact  may  be  used  in  place  of  (2.7)^,  may  be  found  in 
Truesdell  (1969),  p.34,  Eq.  (2.47).  If  the  localization  of  the  Clausius- 
Duhem  inequality  is  instead  carried  out  at  a  point  x^an  the  singular  sur¬ 
face,  the  result  is  the  condition 

fVniO  onS(t),  teftQ.t!]  ,  (3.19) 

which  may  also  be  derived  directly  from  (2.8)3  with  the  help  of  (3.5), 
(3.15)  and  (3.16).  For  the  special  case  of  isothermal,  quasi-static  pro¬ 
cesses  in  thermoelastic  materials,  the  counterpart  of  (3.19)  was  obtained 
by  Knowles  (1979). 

Conversely,  if  the  local  results  (3.18)  and  (3.19)  hold,  then  it  fol¬ 
lows  from  (3. 10) -(3. 12)  that 

rloc(t;D)  +  rcon(t;D)  >  0  ,  r,(t;D)  >  0  ,  (3.20) 

and  hence  from  (3.9)  that  the  Clausius  -  Duhem  inequality  holds.  Also,  it 
is  clear  that  the  entropy  production  rate  T  vanishes  for  every  sub-region  D 
if  and  only  if  equality  holds  in  both  (3.18)  and  (3.19).  If  the  inequality 
in  (3.19)  is  strict  on  Dn8(t),  we  may  conclude  from  (3.9)  and  (3.17)  that 
the  shock  traction  associated  with  the  moving  singular  surface  makes  a  pos¬ 
itive  contribution  to  the  rate  of  entropy  production  T(t;D)  and  thus  repre¬ 
sents  a  dissipative  effect. 
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The  alternate  version  (3.5)  of  the  original  energy  jump  condition 
(2. 8)2  may  be  rewritten  in  still  another  form  with  the  help  of  (3.15)  and 
(3.16): 


P6[[V))  Vn  -  -  f  Vn  -  [  [<^*n]  ]  on  S(t),  te^,^].  (3.21) 

The  original  set  of  jump  conditions  (2.8)  may  be  replaced  by  (2.8)^,  (3.21) 
and  (3.19). 

Note  that  no  constitutive  assumptions  have  yet  been  made.  The  foregoing 
analysis,  however,  would  not  apply  to  the  classical  adiabatic  theory  of 
shock  waves  in  gas  dynamics  (see,  e.g.,  Courant  and  Friedrichs  (1948)), 
since  we  assume  that  the  temperature  is  continuous  across  the  singular 
surface  S(t). 

4.  laalharafli  ataauauLJjCL. a. figntlnuua- 

In  this  section,  we  specialize  the  foregoing  results  to  the  case  of 
isothermal  processes,  by  which  we  mean  those  for  which 

8(x, t)  -  9q  -  constant  for  x€R,  teftQ.tjJ  .  (4.1) 

When  (4.1)  holds,  the  global  form  (2.4)  of  the  first  law  can  be  used  to 
eliminate  the  heat  flux  and  heat  supply  terms  from  the  original  representa¬ 
tion  (2.5)  for  the  entropy  production  rate.  After  doing  so  and  making  use 
of  the  definition  (3.15)  of  the  free  energy  if>,  one  finds  that 
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* 

■ 

*  ’ 

r(t;D)  -  (1/6q) 

an*  v  dA  + 

pb*v  dV  -  d/dt 

p(r/>  +  v*v/2)  dV 

»  « 

3D 

< 

D 

J 

D 

Since  the  contents  of  the  braces  in  (4.2)  represent  the  excess  of  the  rate 
of  mechanical  work  over  the  rate  of  increase  of  the  sum  of  free  and  kinetic 
energies  associated  with  D,  one  infers  from  (4.2)  that,  for  isothermal 
processes,  the  entropy  production  rate  coincides  with  the  rate  of  mechani¬ 
cal  dissipation  per  unit  temperature. 

On  the  other  hand,  with  (4.1)  in  force,  (3.11)  shows  that  the  contrib¬ 
ution  rcon(t;D)  to  T(t;D)  vanishes,  so  that  by  (3.9),  (3.14)  and  (3.17),  we 
also  have 


T(t;D)  -  (1 /60) 


P6  dV  + 

D  S(t)nD 


f  Vn  dA 


Comparison  of  (4.2)  and  (4.3)  yields  the  identity 


(4.3) 


an*v  dA  + 

*+  <*•  a 

3D 


pb«v  dV  + 


(-f)  Vn  dA  - 


s(t)nD 


p6  dV  +  d/dt 
D 


j  P(tf  +  v*v/2)  dV 
D 


for  DCR ,  t€ [  tg ,  t]_ ]  .  (4.4) 
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This  mechanical  balance  motivates  our  choice  of  the  name  shock  driving 
traction  for  f:  Equation  (4.4)  suggests  that  -f(x,t)  be  interpreted  as  the 
component  of  a  vector  traction  acting  normal  to  the  singular  surface  S(t) 
and  exerted  by  the  surface  on  the  body  at  the  point  at  time  t.  Equiva¬ 
lently,  f(x,  t)  may  be  regarded  as  a  normal  traction  applied  to  S(t)  by  the 
body. 


5.  Thermoelastic  materials. 

5.1.  General  processes  in  thermoelastlc  materials.  The  preceding 
discussion  has  made  no  use  of  special  constitutive  relations  for  the  con¬ 
tinuum  under  consideration.  Ve  now  assume  that,  for  the  material  at  hand, 
there  is  a  characterizing  internal  energy  potential  *(F,»j)  such  that 


A 

c(x,t)  -  c  (£(x,  t)  ,»7(x,  t)  )  , 
2,(x,t)  -  />££(£(XIt),»?(x,t)), 
0(x,t)  -  c,(F(x,t)  ,ij(x,t))  ; 


(5.1) 


moreover,  it  is  assumed  that  c^CJF,*)  is  invertible  for  every  tensor  JF  with 
positive  determinant.  We  call  such  a  material  thermoelastic .  To  avoid 
cumbersome  formulas,  we  have  assumed  that  p  and  c(»,«)  are  independent  of 
x,  so  that  the  body  is  homogeneous  in  the  reference  configuration. 


For  a  thermoelastic  material,  (5.1)  imply  that 

pc  -  o*F  +  pSrt  on  R-S(t)  ,  (5.2) 
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so  that  the  internal  dissipation  6  of  (3.13)  vanishes.  It  then  follows  from 
(3.14)  that 


rioc(t;D)  -  0  (5.3) 

for  all  t  and  all  DcR,  whence  (3.9)  becomes 

r(t;D)  -  rcon(t;D)  +  rs(t;D),  (5.4) 

where  rcon  and  Tg  are  given  by  (3.11)  and  (3.17),  respectively.  Thus  in 
every  process  in  a  thermoelastic  material,  entropy  production  is  due 
only  to  heat  conduction  and  the  motion  of  the  singular  surface  S(t). 

By  (5.1)3  and  the  assumed  invertibility  of  £,j(£,  •>,  we  may  write  r)  - 

A  A 

ij(F,9)  and  hence  introduce  the  Helmholtz  free  energy  potential  rf>  for  a 

thermoelastic  material  through 

J(F .»)  -  cCF.nCF.n)  -  *(F, 9)6  .  (5.5) 

Then  by  (3.15),  ^(x.t)  -  ^(F(x, t) , 9 (x, t) ) ,  and  (5.1),  (5.5)  show  that, 
for  any  process , 

o(x,t)  -  p  (F(x,t)  ,9  (x,t))  ,  »?(x,t)  -  -it>g  (F(x,  t)  ,  9  (jt,  t)  )  .  (5.6) 

For  a  thermoelastic  material,  the  shock  driving  traction  f(x,t)  intro- 


duced  for  arbitrary  thermo-mechanical  processes  in  (3.16)  can  be  repre- 

A 

sented  in  terms  of  the  functions  p  and  ^  characteristic  of  the  material  as 
follows : 


f  -  P((J(£,*>n 


<1/2)P 


JF(F,0)  +  $F(F,0) 


[[£]]  . 


(5.7) 


+  +  -  - 

where  F  -  JF(*,t)  and  F  -  F(*,t)  represent  the  limiting  values  of  the 
deformation  gradient  tensor  on  the  positive  and  negative  sides  of  S(t), 
respectively. 


5.2  Isothermal  processes  in  thermoelastic  materials.  For  an  isothermal 
process,  (4.1)  and  (3.11)  show  that  Tcon  -  0  ,  so  that  using  (3.17),  one 
may  further  reduce  (5.4)  to 


r(t;D)  -  rs(t; D) 


(l/*0> 

S(t)nD 


f  vn  dA 


(5.8) 


Thus  for  isothermal  processes  in  thermoelastic  materials,  entropy  produc 
tion  arises  solely  from  the  motion  of  the  surface  of  discontinuity  S(t). 
The  local  consequence  (3.19)  of  the  second  law  at  points  on  the  singular 
surface  of  course  continues  to  hold. 


When  considering  isothermal  processes  in  thermo-elastic  materials,  it 


is  convenient  to  introduce  the 


W  at  the  temperature  6q  by 
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setcing 


W(F)  -  pH f,9q) 


(5.9) 


for  every  tensor  F  with  positive  determinant.  Then  (5.6)^  reads 


o(x,t)  -  Wjr(F(x,  t))  for  xeR-S(t),  te[t0,tj_]. 
The  remaining  mechanical  field  equations  (2. 7)^  2  are 


Div  £  +  pb  -  pv  , 
oFT  -  FaT  , 

#*  P* 


on  R-S(t) , 


(5.10) 


(5.11) 


while  the  mechanical  jump  conditions  are  (2.10)2  and  (2.8)^: 

[[v]]  -  -  C [FJ  J V  ,  [[on]]  +  P[[v]]Vn  -  0  onS(t);  (5.12) 

(5.10)  -  (5.12)  comprise  the  usual  field  equations  and  jump  conditions  of 
elastodynamics .  From  (5.7),  (5.9),  the  shock  driving  traction  f  may  be 
written  in  the  form 


f  -  [[W(£>]]  -  (1/2) 


w£<£)  +  WF<I> 


•[[£]]  . 


(5.13) 


The  inequality  (3.19)  must  also  hold: 


--22-- 


f  Vn  >  0  on  S(t) .  (5.14) 

5.3  Equilibrium  states, and  Isothermal  quasi -static  processes  in  thermo- 
elastic  materials.  Assume  that  the  material  is  thermoelastic,  and  suppose 
that  the  displacement  field  u.  in  (2.1)  is  independent  of  time.  Suppose  fur¬ 
ther  that  the  temperature  6  is  constant,  so  that  (4.1)  holds.  Among  the 
field  equations  (2.7),  the  purely  mechanical  ones  are  now 


Div  a  +  pb  -  0  , 
oFT  -  FerT, 


>■  on  R-S. 


(5.15) 


In  addition,  there  is  the  mechanical  constitutive  requirement  (5.10): 

o  -  Wp(F)  on  R-S.  (5.16) 

The  system  (5.15),  (5.16)  with  F  -  !L  +  Vu  comprises  the  field  equations  of 
elastostatics .  When  this  system  fails  to  be  elliptic,  it  may  possess  weak 
solutions  for  which  £  jumps  across  a  singular  surface  S;  see  Knowles  and 
Sternberg  (1978),  Gurtin  (1983),  Rosakis  (1988).  Such  solutions  must  sat¬ 
isfy  the  purely  mechanical  jump  conditions  (2.10)^  and  (2.8)^  specialized 
to  the  equilibrium  case: 

[knn-O,  [mil-0  on  S.  (5.17) 

***  /v  A** 


J 
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Because  of  the  form  of  the  jump  condition  (5.17)^  in  the  equilibrium 
case,  the  formula  (5.13)  for  the  shock  traction  may  be  simplified  to  give 

♦ 

f-  t[W(F)]]  -  Wjr(F)  •  [  [  F]  ]  ,  (static  case).  (5.18) 

This  representation  of  the  shock  driving  traction  for  equilibrium  deforma¬ 
tions  of  an  elastic  solid  is  equivalent  to  one  derived  by  Yatomi  and  Nish- 
imura  (1983) .  An  alternative  formula  for  f  in  the  present  static  case  is 
readily  shown  to  be 


f  -  n*([P(F)]]n  ,  (static  case).  (5.19) 

where  ii  is  the  unit  normal  to  the  surface  of  discontinuity  S ,  and  P(F) 
is  the  energy -momentum  tensor  introduced  by  Eshelby  (1956,1970,  1975): 

P(F)  -  W<£)1  -  F%£(£)  ;  (5.20) 


here  £  is  the  identity  tensor.  In  the  form  (5 . 19) , (5 . 20)  appropriate  to  the 
equilibrium  case,  the  entity  f  first  appeared  as  the  "force  on  a  defect"  in 
the  work  of  Eshelby  (1956).  This  representation  for  f  is  also  equivalent  to 
that  for  the  force  on  the  interface  between  two  phases  derived  by  Eshelby 
(1970)  and  discussed  by  Rice  (1975).  The  formulas  (5.19),  (5.20)  may  also 
be  found  in  Knowles  (1979). 
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An  isothermal  auasi-static  process  on  a  time  interval  [tQ,t}J  is  a 
one-parameter  family  of  equilibrium  fields  u(*,t)  and  jr(*,t)  with 
smoothness  as  specified  in  Section  2  and  satisfying  (5.15)  -  (5.17)  at  each 
te[tQ,t]J  .  When  at  each  instant  t  the  fields  involve  a  singular  surface 
S(t)  across  which  JF^.t)  jumps,  the  inequality  (5.14)  is  imposed  as  a 
requirement,  f  being  given  by  (5.18)  at  each  t. 

6.  Kinetic  relations. 

Using  the  one - dimens ional ,  isothermal  theory  of  bars  of  non-elliptic 
elastic  material  in  tension  (Abeyaratne  and  Knowles  (1988a, b)),  we  have 
considered  a  bar  lying  along  the  x-axis  in  an  equilibrium  state  in  which 
there  is  a  strain  discontinuity  located  at  an  arbitary  station  x  -  s.  The 
force  F  acting  on  the  bar  is  assumed  to  be  given,  and  the  stress  response 
of  the  material  is  taken  to  be  one  for  which  the  stress  at  first  rises  with 
increasing  strain,  then  declines,  and  finally  rises  again.  For  each  F  in  a 
certain  range,  there  is  a  one  parameter  family,  parameter  s,  of  such 
equilibrium  states  of  physical  interest.  For  each  such  state,  the  overall 
elongation  e  of  the  bar  depends  on  both  F  and  s,  as  do  other  macroscopic 
quantities  such  as  the  total  strain  energy  E  and  the  total  potential  energy 
U  -  E  -  Fe  in  the  bar:  e  -  e(F,s),  E  -  E(F,s),  U  -  U(F,s).  It  turns  out 
that,  in  this  equilibrium  theory,  3U(F,s)/3F  -  -e,  and  3U(F,s)/3s  -  Af  , 
where  A  is  the  cross-sectional  area  of  the  bar,  and  f  is  the  one¬ 
dimensional  counterpart  of  the  static  driving  traction  given  in  (5.18). 

This  suggests  that  the  location  s  of  the  strain  discontinuity  plays  the 
role  of  an  "internal  variable"  whose  "conjugate  force"  is  proportional  to 
the  driving  traction  f.  In  the  references  cited  above,  quasi-static  pro- 
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cesses  are  considered  in  vhich  at  each  instant  the  bar  is  in  one  of  the 
equilibrium  states  just  described.  Admissible  quasi-static  processes  are 
those  that  statisfy  the  one -dimensional  counterpart  fs  >  0  of  (5.14). 

Even  with  admissibility  imposed  as  a  requirement,  specifying  the  force  his¬ 
tory  F (r),  0  S  r  :S  t,  in  a  quasi-static  motion  fails  to  determine  the 
current  value  of  the  elongation  e(t),  since  the  shock  location  (or  internal 
variable)  s(t),  must  be  specified  as  well.  This  suggests  a  constitutive 
deficiency.  In  internal  variable  theories  designed  to  model  microstructural 
effects  on  inelastic  macro structural  behavior,  a  so-called  "kinetic  rela¬ 
tion"  giving  the  time  rate  of  change  of  the  internal,  or  microstructural, 
variable  as  a  function  of  the  associated  "thermodynamic  force"  is  commonly 
added  as  a  part  of  the  constitutive  description;  see,  for  example,  Rice 
(1970,  1971,  1975).  In  their  one -dimensional  setting,  Abeyaratne  and 
Knowles  ( 1988a, b)  adopt  this  point  of  view  and  hence  require  that  the  shock 
location  s(t)  and  the  associated  driving  traction  f(t)  be  related  by  a 
kinetic  law  of  the  form  s(t)  -  V(f(t)),  where  the  kinetic  response  function 
V  is  determined  by  the  material.  Such  a  kinetic  relation  then  couples  the 
time -evolution  of  the  location  of  the  surface  of  strain  discontinuity  to 
the  local  strains  on  either  side  of  the  jump. 

The  isothermal  equilibrium  problem  of  the  twisting  of  an  infinite 
medium  containing  a  circular  hole  considered  by  Abeyaratne  and  Knowles 
(1987a, b)  involves  cylindrical  geometry  in  plane  strain  and  is  essentially 
one -dimensional  because  of  axial  symmetry.  The  surface  of  discontinuity  is 
now  circular,  and  its  radius  plays  the  role  of  the  location  s  of  the  strain 
jump  in  the  tensile  bar  problem.  The  formalism  appropriate  to  internal 
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variables  emerges  in  this  problem  as  well,  and  an  additional  constitutive 
postulate  in  the  form  of  a  kinetic  relation  may  be  proposed  here  on  similar 
grounds.  It  would  now  take  the  form 

Vn(x,t)  -  V(f(x,t)),  x€S(t).  (6.1) 

where  x_ is  the  position  vector  of  a  typical  point  on  the  cylindrical  sur¬ 
face  of  discontinuity  S(t),  f(x,t)  is  given  at  each  t  by  (5.18)  specialized 
to  plane  strain  for  incompressible  materials,  and  Vn(x,t)  is  the  component 
of  velocity  of  the  singular  surface  normal  to  itself,  i.e.  in  the  radial 
direction.  Again,  V(f)  is  a  function  determined  by  the  material.  To  be  con¬ 
sistent  with  the  admissibility  requirement  (5.14),  V  must  be  such  that 

V(f)f  £  0  (6.2) 

for  all  possible  values  of  f. 

The  consequences  of  imposing  a  kinetic  relation  in  the  one -dimensional 
setting  of  the  bar  are  discussed  in  detail  in  Abeyaratne  and  Knowles 
(1988a, b)  and  will  not  be  repeated  here.  A  similar  discussion  could  be  pro¬ 
vided  for  the  twist  problem.  While  in  both  cases  the  motivation  for  the 
imposition  of  the  kinetic  relation  as  a  requirement  is  taken  from  internal 
variable  theories  of  inelastic  behavior,  the  possibility  of  imposing  it 
comes  about  because  of  the  lack  of  uniqueness  of  weak  solutions  to  the 
relevant  equilibrium  problem.  Indeed,  for  a  given  force  F  in  the  bar  prob¬ 
lem,  for  example,  solutions  with  a  single  "equilibrium  shock"  at  x  -  s  fail 
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to  be  unique  precisely  Co  the  extent  that  their  totality  comprises  a  one- 
parameter  family  (parameter  s) .  This  is  exactly  the  extent  of  lack-of- 
uniqueness  needed  to  make  room  for  the  kinetic  law  as  an  additional 
requirement.  Thus  the  equilibrium  problem  seems  to  be  under -determined  to 
precisely  the  degree  necessary  to  accommodate  a  kinetic  relation  applicable 
to  isothermal  ouasl-static  processes  in  the  bar.  Entirely  similar  remarks 
apply  to  the  axially  symmetric  twist  problem. 

The  general  proposal  of  a  constitutive  relationship  between  the  driving 
traction  f(jc, t)  and  the  normal  velocity  Vn(x, t)  of  the  material  singular 
surface  S(t)  during  quasi-static  processes  in  a  thermoelastic  material  may 
be  examined  from  the  perspective  of  irreversible  thermodynamics.  In  the 
representation  for  the  entropy  production  rate  r  provided  by  (5.4),  (3.11) 
and  (3.17),  one  could  view  the  ratios  (Grad  9)/0 ^  and  f/$  as  thermodynamic 
"affinities"  and  the  terms  q  and  V_  as  the  corresponding  "fluxes";  see 
Chapter  14  of  Callen  (1985),  Chapter  14  of  Kestin  (1968),  and  Lecture  7  of 
Truesdell  (1969)  for  discussions  of  these  notions.  In  the  theory  of  irrev- 
erslble  processes,  it  is  customary  to  postulate  a  constitutive  relationship 
in  which  the  present  value  of  each  flux  is  a  function  of  the  present  value 
of  the  affinities  and  perhaps  of  their  past  histories  as  well.  In  our  set¬ 
ting,  one  might  -  as  a  simplest  case  -  postulate  a  relationship  between  the 
present  values  of  each  flux  and  its  corresponding  affinity,  as  in  the 
theory  of  "purely  resistive"  thermodynamical  systems  (Callen  (1985),  Chap¬ 
ter  14) .  This  would  give 

q(x,t)  -  Q*(Grad  0(x, t)/02(x,t)) ,  xeR-S(t), 


(6.3) 
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and 

Vn(£,t)  -  V*(f(x,t)/*(x,t)),  xeS(t)  ,  (6.4) 

where  and  V*  are  functions  determined  by  the  material;  the  constitutive 
statements  (6.3)  and  (6.4)  represent  a  heat  conduction  law  and  a  kinetic 
relation,  respectively.  The  consequence  (5.14)  of  the  Claus ius-Duhem 
inequality  would  require  that 

V*(f/rf)f/*  2:  0  (6.5) 

for  all  possible  values  of  f/B .  For  isothermal  processes  of  the  type 
discussed  in  the  preceding  paragraphs  of  this  section,  (6.4)  reduces  to 
(6.1)  with  V(f)  -  V*(f /Bq),  and  (6.5)  reduces  to  (6.2).  We  observe  that,  if 
the  kinetic  response  function  V*  (or,  equivalently,  V)  is  continuous,  then 
(6.5)  requires  that 

V*(0)  -  V(0)  -  0.  (6.6) 

If  the  body  is  at  a  uniform  temperature  and  in  a  state  of  mechanical 
equilibrium  involving  an  equilibrium  shock  S  at  all  points  of  which  the 
Maxwell  condition  f  —  0  holds,  we  say  the  body  is  in  a  Maxwell  state.  If  V* 
is  smooth  and  departures  from  a  Maxwell  state  are  slight,  so  that  f/6  is 
small,  one  might  replace  (6.4)  by  its  linearization.  (Linearized  kinetic 
relations  are  often  used  in  irreversible  thermodynamics  to  describe  pro- 
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cesses  thaC  are  "close  to  thermodynamic  equilibrium".)  In  view  of  (6.6), 
this  would  yield 

Vn(x,t)  -  »/f(x,t)/0(x,t)  ,  xeS(t),  (6.7) 

where  v  -  V'(0)  Is  a  material  constant;  by  (6.2),  v  is  necessarily  posi¬ 
tive.  A  linear  relation  of  the  type  (6.7)  between  the  driving  traction  and 
the  velocity  of  the  singular  surface  was  shown  by  Abeyaratne  and  Knowles 
(1988b)  to  lead  to  a  conventional  type  of  viscoelastic  macroscopic  response 
in  the  one -dimensional  theory  of  quasi-static  processes  in  tensile  bars 
composed  of  a  particular  non- elliptic  elastic  material. 

Note  that  the  singular  surface  S(t)  lies  in  the  region  occupied  by  the 
body  in  the  reference  configuration;  thus  the  driving  traction  f  might  more 
precisely  be  called  the  nominal  driving  traction.  The  kinetic  relation 
(6.4)  is  thus  a  "Lagrangian" ,  or  "material",  assertion.  It  is  readily  shown 
that  it  conforms  to  the  principle  of  material  frame  indifference. 

For  isothermal,  quasi-static  processes  in  a  thermoelastic  material,  the 
state  of  the  body  at  each  instant  is  one  of  mechanical  equilibrium,  and  an 
equilibrium  shock  cannot  occur  unless  the  corresponding  elastic  potential 
is  non-elliptic .  The  equilibrium  shock  may  be  thought  of  as  a  phase  bound¬ 
ary,  and  the  quasi-static  process  may  be  viewed  as  one  in  which  particles 
of  the  body  are  being  transformed  from  one  phase  to  another.  The  normal 
velocity  Vn  is  clearly  a  measure  of  the  rate  at  which  this  takes  place,  and 
the  relation  (6.4)  may  thus  be  regarded  as  controlling  the  kinetics  of  the 
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phase  transformation;  see,  for  example.  Section  3.1  of  Fine  (1964)  or  Sec¬ 
tion  1.9  of  Porter  and  Easterling  (1981)  for  related  discussions. 

If  the  kinetic  relation  is  to  be  constitutive,  it  must  apply  not  only 
during  quasi-static  processes,  but  when  inertial  effects  are  included  as 
well.  For  dynamical  processes,  it  is  to  be  expected  that  moving  surfaces  of 
strain  discontinuity  S(t)  -  counterparts  of  ordinary  shock  waves,  for 
example  -  can  occur  which  are  not  phase  boundaries ,  and  for  which  a 
supplementary  relation  such  as  (6.4)  is  thus  expected  to  be  inappropriate 
on  physical  grounds.  From  the  mathematical  point  of  view,  it  is  the  lack 
of  uniqueness  of  solutions  of  equilibrium  boundary  value  problems  for 
non-elliptic  elastic  materials  that  allows  the  prescription  of  a  supplemen¬ 
tary  kinetic  relation  for  quasi-static  processes.  Whether  there  is  a  pre¬ 
cisely  analogous  lack  of  uniqueness  for  the  boundary- initial-value  problems 
arising  in  the  dynamics  of  non-elliptic  elastic  materials,  and  how  in 
general  to  distinguish  phase  boundaries  from  ordinary  shock  waves,  become 
questions  of  some  importance.  In  the  mathematical  study  of  systems  of  con¬ 
servation  laws  in  one  space  dimension  (see,  for  example,  Dafermos 
(1983,1984),  Lax  (1973)),  it  is  known  that  the  solution  to  the  initial- 
value  problem  subject  to  an  "entropy  inequality"  such  as  (2.8)3  (or,  equi¬ 
valently,  (3.19))  is  unique,  provided  that  the  curvature  of  the  underlying 
stress-strain  relation  is  always  of  one  sign.  In  the  absence  of  both  con¬ 
vexity  and  concavity,  the  entropy  inequality  is  not  strong  enough  to  secure 
uniqueness.  The  fact  that  the  entropy  inequality  must  be  supplemented  in 
order  to  assure  uniqueness  for  the  initial  value  problem  for  such  materials 
reflects  the  need  for  an  additional  requirement  of  a  constitutive  nature, 


--31-- 


as  distinguished  from  a  fundamental  thermo -mechanical  principle.  A  kinetic 
relation  such  as  (6.4),  in  which  f  is  now  given  by  the  representation  (5.7) 
pertinent  to  dynamical  processes  and  the  material  function  V*  is  subject  to 
the  requirement  (6.5)  imposed  by  the  entropy  inequality,  might  provide  an 
appropriate  supplementary  constitutive  requirement.  Additional  restric¬ 
tions  strong  enough  to  guarantee  uniqueness  in  the  initial  value  problem  in 
one  space  dimension  for  materials  whose  stress- strain  relation  is  neither 
convex  nor  concave  have  been  investigated  from  viewpoints  different  from 
the  present  one  by  a  number  of  authors:  see,  for  example,  Hatori  (1986), 
James  (1980),  Oleinik  (1959),  Shearer  (1985)  and  Slemrod  (1983). 
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